Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis.

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

Shoring up the base: the development and regulation of cortical sclerenchyma in grass nodal roots.

Plants depend on the combined action of a shoot-root-soil system to maintain their anchorage to the soil. Mechanical failure of any component of this system results in lodging, a permanent and irreversible inability to maintain vertical orientation. Models of anchorage in grass crops identify the compressive strength of roots near the soil surface as key determinant of resistance to lodging. Indeed, studies of disparate grasses report a ring of thickened, sclerenchyma cells surrounding the root cortex, present only at the base of nodal roots. Here, in the investigation of the development and regulation of this agronomically important trait, we show that development of these cells is uncoupled from the maturation of other secondary cell wall-fortified cells, and that cortical sclerenchyma wall thickening is stimulated by mechanical forces transduced from the shoot to the root. We also show that exogenous application of gibberellic acid stimulates thickening of lignified cell types in the root, including cortical sclerenchyma, but is not sufficient to establish sclerenchyma identity in cortex cells. Leveraging the ability to manipulate cortex development via mechanical stimulus, we show that cortical sclerenchyma development alters root mechanical properties and improves resistance to lodging. We describe transcriptome changes associated with cortical sclerenchyma development under both ambient and mechanically stimulated conditions and identify SECONDARY WALL NAC7 as a putative regulator of mechanically responsive cortex cell wall development at the root base.

Functional analysis of Salix purpurea genes support roles for ARR17 and GATA15 as master regulators of sex determination.

The Salicaceae family is of growing interest in the study of dioecy in plants because the sex determination region (SDR) has been shown to be highly dynamic, with differing locations and heterogametic systems between species. Without the ability to transform and regenerate Salix in tissue culture, previous studies investigating the mechanisms regulating sex in the genus Salix have been limited to genome resequencing and differential gene expression, which are mostly descriptive in nature, and functional validation of candidate sex determination genes has not yet been conducted. Here, we used Arabidopsis to functionally characterize a suite of previously identified candidate genes involved in sex determination and sex dimorphism in the bioenergy shrub willow Salix purpurea. Six candidate master regulator genes for sex determination were heterologously expressed in Arabidopsis, followed by floral proteome analysis. In addition, 11 transcription factors with predicted roles in mediating sex dimorphism downstream of the SDR were tested using DAP-Seq in both male and female S. purpurea DNA. The results of this study provide further evidence to support models for the roles of ARR17 and GATA15 as master regulator genes of sex determination in S. purpurea, contributing to a regulatory system that is notably different from that of its sister genus Populus. Evidence was also obtained for the roles of two transcription factors, an AP2/ERF family gene and a homeodomain-like transcription factor, in downstream regulation of sex dimorphism.

A qualitative meta-analysis exploring client-reported outcomes of couple therapy.

The quantitative reviews of the outcome research on couple therapy show that this type of therapy can produce positive outcomes for couples and improve relationship satisfaction. There is now also a number of qualitative studies in which clients report in their own words on the outcomes of couple therapy. This study aimed to meta-analyze the client-reported outcomes of couple therapy generated in the studies using qualitative methods. A sample of 15 primary studies examining clients' reported outcomes of couple therapy was identified through an extensive literature search. Relevant qualitative data on the client-reported outcomes were extracted into a single data set. The data was then analyzed using a descriptive-interpretive qualitative meta-analytic approach. Similar outcomes were grouped into metacategories. The metacategories were then organized into several clusters of the client-reported outcomes of couple therapy. The meta-analysis yielded 25 metacategories which were clustered into seven main clusters, (a) seeing things differently; (b) changed behavior within the relationship; (c) improved experience in the relationship; (d) improved communication quality; (e) improvement in relationship functioning; (f) improved individual functioning, and (g) difficult outcomes of therapy. Clients reported numerous constructive (e.g., new understanding of the couple's interactional functioning, improvement in the conflict management, new positive ways of relating and connecting, letting go of expectations imposed on the partner or changes within the self that may be one's contribution to the relationship), and some difficult, outcomes of engaging in couple therapy (the clarity on the decision to separate). (PsycInfo Database Record (c) 2023 APA, all rights reserved).

Speech patterns in responses to questions asked by an intelligent virtual agent can help to distinguish between people with early stage neurodegenerative disorders and healthy controls.

Previous research has provided strong evidence that speech patterns can help to distinguish between people with early stage neurodegenerative disorders (ND) and healthy controls. This study examined speech patterns in responses to questions asked by an intelligent virtual agent (IVA): a talking head on a computer which asks pre-recorded questions. The study investigated whether measures of response length, speech rate and pausing in responses to questions asked by an IVA help to distinguish between healthy control participants and people diagnosed with Mild Cognitive Impairment (MCI) or Alzheimer's disease (AD). The study also considered whether those measures can further help to distinguish between people with MCI, people with AD, and healthy control participants (HC). There were 38 people with ND (31 people with MCI, 7 people with AD) and 26 HC. All interactions took place in English. People with MCI spoke fewer words compared to HC, and people with AD and people with MCI spoke for less time than HC. People with AD spoke at a slower rate than people with MCI and HC. There were significant differences across all three groups for the proportion of time spent pausing and the average pause duration: silent pauses make up the greatest proportion of responses from people with AD, who also have the longest average silent pause duration, followed by people with MCI then HC. Therefore, the study demonstrates the potential of an IVA as a method for collecting data showing patterns which can help to distinguish between diagnostic groups.

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning.

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies-Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR's ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods. MetaBAT-LR is open-source and is available at https://bitbucket.org/project-metabat/metabat-lr.

IgG4 Disease-Related Ataxia.

We describe a male patient presenting with cerebellar ataxia and behavioural frontotemporal dementia in whom imaging showed cerebellar atrophy. He had significantly low N-acetyl aspartate to creatine (NAA/Cr) area ratio on MR spectroscopy of the cerebellum, primarily affecting the vermis. CT body scan showed extensive abnormal tissue within the mesentery, the retroperitoneum and perinephric areas. PET-CT showed increased tracer uptake within the wall of the aorta suggestive of an aortitis and within the perinephric tissue bilaterally. Biopsy of the perinephric tissue confirmed IgG4 disease. Treatment with steroids and mycophenolate improved his clinical state, but he developed symptoms attributed to pericardiac effusion that necessitated treatment initially with drainage and subsequently with pericardial window. After a course of rituximab, he had an episode of sepsis that did not respond to appropriate treatment and died as a result. Both the imaging findings and neurological presentation with cerebellar ataxia and behavioural frontotemporal dementia are novel in the context of IgG4 disease.

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

High-throughput identification of novel heat tolerance genes via genome-wide pooled mutant screens in the model green alga Chlamydomonas reinhardtii.

Different high temperatures adversely affect crop and algal yields with various responses in photosynthetic cells. The list of genes required for thermotolerance remains elusive. Additionally, it is unclear how carbon source availability affects heat responses in plants and algae. We utilized the insertional, indexed, genome-saturating mutant library of the unicellular, eukaryotic green alga Chlamydomonas reinhardtii to perform genome-wide, quantitative, pooled screens under moderate (35°C) or acute (40°C) high temperatures with or without organic carbon sources. We identified heat-sensitive mutants based on quantitative growth rates and identified putative heat tolerance genes (HTGs). By triangulating HTGs with heat-induced transcripts or proteins in wildtype cultures and MapMan functional annotations, we presented a high/medium-confidence list of 933 Chlamydomonas genes with putative roles in heat tolerance. Triangulated HTGs include those with known thermotolerance roles and novel genes with little or no functional annotation. About 50% of these high-confidence HTGs in Chlamydomonas have orthologs in green lineage organisms, including crop species. Arabidopsis thaliana mutants deficient in the ortholog of a high-confidence Chlamydomonas HTG were also heat sensitive. This work expands our knowledge of heat responses in photosynthetic cells and provides engineering targets to improve thermotolerance in algae and crops.

Genome of Paspalum vaginatum and the role of trehalose mediated autophagy in increasing maize biomass.

A number of crop wild relatives can tolerate extreme stress to a degree outside the range observed in their domesticated relatives. However, it is unclear whether or how the molecular mechanisms employed by these species can be translated to domesticated crops. Paspalum (Paspalum vaginatum) is a self-incompatible and multiply stress-tolerant wild relative of maize and sorghum. Here, we describe the sequencing and pseudomolecule level assembly of a vegetatively propagated accession of P. vaginatum. Phylogenetic analysis based on 6,151 single-copy syntenic orthologues conserved in 6 related grass species places paspalum as an outgroup of the maize-sorghum clade. In parallel metabolic experiments, paspalum, but neither maize nor sorghum, exhibits a significant increase in trehalose when grown under nutrient-deficit conditions. Inducing trehalose accumulation in maize, imitating the metabolic phenotype of paspalum, results in autophagy dependent increases in biomass accumulation.

The landscape of Chlamydomonas histone H3 lysine 4 methylation reveals both constant features and dynamic changes during the diurnal cycle.

Chromatin modifications are epigenetic regulatory features with major roles in various cellular events, yet they remain understudied in algae. We interrogated the genome-wide distribution pattern of mono- and trimethylated histone H3 lysine 4 (H3K4) using chromatin-immunoprecipitation followed by deep-sequencing (ChIP-seq) during key phases of the Chlamydomonas cell cycle: early G1 phase, Zeitgeber Time 1 (ZT1), when cells initiate biomass accumulation, S/M phase (ZT13) when cells are replicating DNA and undergoing mitosis, and late G0 phase (ZT23) when they are quiescent. Tri-methylated H3K4 was predominantly enriched at transcription start sites of the majority of protein coding genes (85%). The likelihood of a gene being marked by H3K4me3 correlated with it being transcribed at some point during the life cycle but not necessarily by continuous active transcription, as exemplified by early zygotic genes, which may remain transcriptionally dormant for thousands of generations between sexual cycles. The exceptions to this rule were around 120 loci, some of which encode non-poly-adenylated transcripts, such as small nuclear RNAs and replication-dependent histones that had H3K4me3 peaks only when they were being transcribed. Mono-methylated H3K4 was the default state for the vast majority of histones that were bound outside of transcription start sites and terminator regions of genes. A small fraction of the genome that was depleted of any H3 lysine 4 methylation was enriched for DNA cytosine methylation and the genes within these DNA methylation islands were poorly expressed. Besides marking protein coding genes, H3K4me3 ChIP-seq data served also as a annotation tool for validation of hundreds of long non-coding RNA genes.

Sodalis ligni Strain 159R Isolated from an Anaerobic Lignin-Degrading Consortium.

Novel bacterial isolates with the capabilities of lignin depolymerization, catabolism, or both, could be pertinent to lignocellulosic biofuel applications. In this study, we aimed to identify anaerobic bacteria that could address the economic challenges faced with microbial-mediated biotechnologies, such as the need for aeration and mixing. Using a consortium seeded from temperate forest soil and enriched under anoxic conditions with organosolv lignin as the sole carbon source, we successfully isolated a novel bacterium, designated 159R. Based on the 16S rRNA gene, the isolate belongs to the genus Sodalis in the family Bruguierivoracaceae. Whole-genome sequencing revealed a genome size of 6.38 Mbp and a GC content of 55 mol%. To resolve the phylogenetic position of 159R, its phylogeny was reconstructed using (i) 16S rRNA genes of its closest relatives, (ii) multilocus sequence analysis (MLSA) of 100 genes, (iii) 49 clusters of orthologous groups (COG) domains, and (iv) 400 conserved proteins. Isolate 159R was closely related to the deadwood associated Sodalis guild rather than the tsetse fly and other insect endosymbiont guilds. Estimated genome-sequence-based digital DNA-DNA hybridization (dDDH), genome percentage of conserved proteins (POCP), and an alignment analysis between 159R and the Sodalis clade species further supported that isolate 159R was part of the Sodalis genus and a strain of Sodalis ligni. We proposed the name Sodalis ligni str. 159R (=DSM 110549 = ATCC TSD-177). IMPORTANCE Currently, in the paper industry, paper mill pulping relies on unsustainable and costly processes to remove lignin from lignocellulosic material. A greener approach is biopulping, which uses microbes and their enzymes to break down lignin. However, there are limitations to biopulping that prevent it from outcompeting other pulping processes, such as requiring constant aeration and mixing. Anaerobic bacteria are a promising alternative source for consolidated depolymerization of lignin and its conversion to valuable by-products. We presented Sodalis ligni str. 159R and its characteristics as another example of potential mechanisms that can be developed for lignocellulosic applications.

An optimized ChIP-Seq framework for profiling histone modifications in Chromochloris zofingiensis.

The eukaryotic green alga Chromochloris zofingiensis is a reference organism for studying carbon partitioning and a promising candidate for the production of biofuel precursors. Recent transcriptome profiling transformed our understanding of its biology and generally algal biology, but epigenetic regulation remains understudied and represents a fundamental gap in our understanding of algal gene expression. Chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) is a powerful tool for the discovery of such mechanisms, by identifying genome-wide histone modification patterns and transcription factor-binding sites alike. Here, we established a ChIP-Seq framework for Chr. zofingiensis yielding over 20 million high-quality reads per sample. The most critical steps in a ChIP experiment were optimized, including DNA shearing to obtain an average DNA fragment size of 250 bp and assessment of the recommended formaldehyde concentration for optimal DNA-protein cross-linking. We used this ChIP-Seq framework to generate a genome-wide map of the H3K4me3 distribution pattern and to integrate these data with matching RNA-Seq data. In line with observations from other organisms, H3K4me3 marks predominantly transcription start sites of genes. Our H3K4me3 ChIP-Seq data will pave the way for improved genome structural annotation in the emerging reference alga Chr. zofingiensis.

The role of ATXR6 expression in modulating genome stability and transposable element repression in Arabidopsis.

ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) AND ATXR6 are required for the deposition of H3K27me1 and for maintaining genomic stability in Arabidopsis Reduction of ATXR5/6 activity results in activation of DNA damage response genes, along with tissue-specific derepression of transposable elements (TEs), chromocenter decompaction, and genomic instability characterized by accumulation of excess DNA from heterochromatin. How loss of ATXR5/6 and H3K27me1 leads to these phenotypes remains unclear. Here we provide extensive characterization of the atxr5/6 hypomorphic mutant by comprehensively examining gene expression and epigenetic changes in the mutant. We found that the tissue-specific phenotypes of TE derepression and excessive DNA in this atxr5/6 mutant correlated with residual ATXR6 expression from the hypomorphic ATXR6 allele. However, up-regulation of DNA damage genes occurred regardless of ATXR6 levels and thus appears to be a separable process. We also isolated an atxr6-null allele which showed that ATXR5 and ATXR6 are required for female germline development. Finally, we characterize three previously reported suppressors of the hypomorphic atxr5/6 mutant and show that these rescue atxr5/6 via distinct mechanisms, two of which involve increasing H3K27me1 levels.

Persistence and plasticity in bacterial gene regulation.

Organisms orchestrate cellular functions through transcription factor (TF) interactions with their target genes, although these regulatory relationships are largely unknown in most species. Here we report a high-throughput approach for characterizing TF-target gene interactions across species and its application to 354 TFs across 48 bacteria, generating 17,000 genome-wide binding maps. This dataset revealed themes of ancient conservation and rapid evolution of regulatory modules. We observed rewiring, where the TF sensing and regulatory role is maintained while the arrangement and identity of target genes diverges, in some cases encoding entirely new functions. We further integrated phenotypic information to define new functional regulatory modules and pathways. Finally, we identified 242 new TF DNA binding motifs, including a 70% increase of known Escherichia coli motifs and the first annotation in Pseudomonas simiae, revealing deep conservation in bacterial promoter architecture. Our method provides a versatile tool for functional characterization of genetic pathways in prokaryotes and eukaryotes.

Long-read metagenomics of soil communities reveals phylum-specific secondary metabolite dynamics.

Microbial biosynthetic gene clusters (BGCs) encoding secondary metabolites are thought to impact a plethora of biologically mediated environmental processes, yet their discovery and functional characterization in natural microbiomes remains challenging. Here we describe deep long-read sequencing and assembly of metagenomes from biological soil crusts, a group of soil communities that are rich in BGCs. Taking advantage of the unusually long assemblies produced by this approach, we recovered nearly 3,000 BGCs for analysis, including 712 full-length BGCs. Functional exploration through metatranscriptome analysis of a 3-day wetting experiment uncovered phylum-specific BGC expression upon activation from dormancy, elucidating distinct roles and complex phylogenetic and temporal dynamics in wetting processes. For example, a pronounced increase in BGC transcription occurs at night primarily in cyanobacteria, implicating BGCs in nutrient scavenging roles and niche competition. Taken together, our results demonstrate that long-read metagenomic sequencing combined with metatranscriptomic analysis provides a direct view into the functional dynamics of BGCs in environmental processes and suggests a central role of secondary metabolites in maintaining phylogenetically conserved niches within biocrusts.

Aspects of the Neurospora crassa Sulfur Starvation Response Are Revealed by Transcriptional Profiling and DNA Affinity Purification Sequencing.

Accurate nutrient sensing is important for rapid fungal growth and exploitation of available resources. Sulfur is an important nutrient source found in a number of biological macromolecules, including proteins and lipids. The model filamentous fungus Neurospora crassa is capable of utilizing sulfur found in a variety of sources from amino acids to sulfate. During sulfur starvation, the transcription factor CYS-3 is responsible for upregulation of genes involved in sulfur uptake and assimilation. Using a combination of RNA sequencing and DNA affinity purification sequencing, we performed a global survey of the N. crassa sulfur starvation response and the role of CYS-3 in regulating sulfur-responsive genes. The CYS-3 transcription factor bound the promoters and regulated genes involved in sulfur metabolism. Additionally, CYS-3 directly activated the expression of a number of uncharacterized transporter genes, suggesting that regulation of sulfur import is an important aspect of regulation by CYS-3. CYS-3 also directly regulated the expression of genes involved in mitochondrial electron transfer. During sulfur starvation, genes involved in nitrogen metabolism, such as amino acid and nucleic acid metabolic pathways, along with genes encoding proteases and nucleases that are necessary for scavenging nitrogen, were activated. Sulfur starvation also caused changes in the expression of genes involved in carbohydrate metabolism, such as those encoding glycosyl hydrolases. Thus, our data suggest a connection between sulfur metabolism and other aspects of cellular metabolism. IMPORTANCE Identification of nutrients present in the environment is a challenge common to all organisms. Sulfur is an important nutrient source found in proteins, lipids, and electron carriers that are required for the survival of filamentous fungi such as Neurospora crassa. Here, we transcriptionally profiled the response of N. crassa to characterize the global response to sulfur starvation. We also used DNA affinity purification sequencing to identify the direct downstream targets of the transcription factor responsible for regulating genes involved in sulfur uptake and assimilation. Along with genes involved in sulfur metabolism, this transcription factor regulated a number of uncharacterized transporter genes and genes involved in mitochondrial electron transfer. Our data also suggest a connection between sulfur, nitrogen, and carbon metabolism, indicating that the regulation of a number of metabolic pathways is intertwined.

Plant single-cell solutions for energy and the environment.

Progress in sequencing, microfluidics, and analysis strategies has revolutionized the granularity at which multicellular organisms can be studied. In particular, single-cell transcriptomics has led to fundamental new insights into animal biology, such as the discovery of new cell types and cell type-specific disease processes. However, the application of single-cell approaches to plants, fungi, algae, or bacteria (environmental organisms) has been far more limited, largely due to the challenges posed by polysaccharide walls surrounding these species' cells. In this perspective, we discuss opportunities afforded by single-cell technologies for energy and environmental science and grand challenges that must be tackled to apply these approaches to plants, fungi and algae. We highlight the need to develop better and more comprehensive single-cell technologies, analysis and visualization tools, and tissue preparation methods. We advocate for the creation of a centralized, open-access database to house plant single-cell data. Finally, we consider how such efforts should balance the need for deep characterization of select model species while still capturing the diversity in the plant kingdom. Investments into the development of methods, their application to relevant species, and the creation of resources to support data dissemination will enable groundbreaking insights to propel energy and environmental science forward.

DOE JGI Metagenome Workflow.

The DOE Joint Genome Institute (JGI) Metagenome Workflow performs metagenome data processing, including assembly; structural, functional, and taxonomic annotation; and binning of metagenomic data sets that are subsequently included into the Integrated Microbial Genomes and Microbiomes (IMG/M) (I.-M. A. Chen, K. Chu, K. Palaniappan, A. Ratner, et al., Nucleic Acids Res, 49:D751-D763, 2021, https://doi.org/10.1093/nar/gkaa939) comparative analysis system and provided for download via the JGI data portal (https://genome.jgi.doe.gov/portal/). This workflow scales to run on thousands of metagenome samples per year, which can vary by the complexity of microbial communities and sequencing depth. Here, we describe the different tools, databases, and parameters used at different steps of the workflow to help with the interpretation of metagenome data available in IMG and to enable researchers to apply this workflow to their own data. We use 20 publicly available sediment metagenomes to illustrate the computing requirements for the different steps and highlight the typical results of data processing. The workflow modules for read filtering and metagenome assembly are available as a workflow description language (WDL) file (https://code.jgi.doe.gov/BFoster/jgi_meta_wdl). The workflow modules for annotation and binning are provided as a service to the user community at https://img.jgi.doe.gov/submit and require filling out the project and associated metadata descriptions in the Genomes OnLine Database (GOLD) (S. Mukherjee, D. Stamatis, J. Bertsch, G. Ovchinnikova, et al., Nucleic Acids Res, 49:D723-D733, 2021, https://doi.org/10.1093/nar/gkaa983).IMPORTANCE The DOE JGI Metagenome Workflow is designed for processing metagenomic data sets starting from Illumina fastq files. It performs data preprocessing, error correction, assembly, structural and functional annotation, and binning. The results of processing are provided in several standard formats, such as fasta and gff, and can be used for subsequent integration into the Integrated Microbial Genomes and Microbiomes (IMG/M) system where they can be compared to a comprehensive set of publicly available metagenomes. As of 30 July 2020, 7,155 JGI metagenomes have been processed by the DOE JGI Metagenome Workflow. Here, we present a metagenome workflow developed at the JGI that generates rich data in standard formats and has been optimized for downstream analyses ranging from assessment of the functional and taxonomic composition of microbial communities to genome-resolved metagenomics and the identification and characterization of novel taxa. This workflow is currently being used to analyze thousands of metagenomic data sets in a consistent and standardized manner.